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(Dated: July 2, 2007) 

Abstract 

Using 281 pb" 1 of data recorded by the CLEO-c detector in e + e~ collisions at the ^(3770), 
corresponding to 0.78 million D + D~ pairs, we investigate the substructure of the decay D + — > 
7r~7r + 7r + using the Dalitz plot technique. We find that our data are consistent with the following 
intermediate states: /9(770)vr+, / 2 (1270)7r+, /o(1370)vr+, / (1500)vr+, /o(980)vr + , and air + . We 
confirm large S wave contributions at low irir mass. We set upper limits on contributions of other 
possible intermediate states. We consider three models of the tttt S wave and find that all of them 
adequately describe our data. 
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I. INTRODUCTION 



The study of charmed meson hadronic decays illuminates light meson spectroscopy. Many 
of these decays proceed via quasi two-body modes and are subsequently observed as three 
or more stable particles. In this work our goal is to describe the two-body resonances that 
contribute to the observed three-body D + — > tt~tt + tt + decay. Study of a given state can 
shed light on different production mechanisms. 

We present here a study of charged D decay to three charged pions carried out with the 
CLEO detector. This mode has been studied previously by E687 [l], E691 0], E791 0], and 
FOCUS The analyses from E791 and FOCUS have roughly the same data size as the 
one described here, while the E687 and E691 analyses used about an order of magnitude 
smaller samples and are not discussed further. 

E791 uses the isobar technique, where each resonant contribution to the Dalitz plot [B| 
is modeled as a Breit-Wigner amplitude with a complex phase. This works well for narrow, 
well separated resonances, but when the resonances are wide and start to overlap, solutions 
become ambiguous, and unitarity is violated. In contrast, FOCUS uses the K-matrix ap- 
proach, which gives a description of S wave tttt resonances treating the a (also known as 
/o(600)) and Al(980) contributions in a unified way. While this approach is a step forward, 
some authors jfjj, 0] have claimed that the exact formalism used by FOCUS violates chiral 
constraints, and might therefore lead to unphysical behavior at low tttt mass, where the S 
wave is most prominent. Despite the difference in approach the two techniques give a good 
description of the observed Dalitz plots and agree about the overall contributions of the res- 
onances, as is shown in Table [B Both experiments see that about half of the fit fraction for 
this decay is explained by a low tt + tt~ mass S wave. We have in hand a comparable sample 
of D + — > tt~tt + tt + decays (inclusion of the charge-conjugate mode is always implicit); we 
can thus check this somewhat surprising result in a significantly different environment. 

E791 and FOCUS are fixed target experiments where D mesons are produced within 
a momentum range of 10-100 GeV/c. In our experiment D + mesons are produced in the 
process e + e~ — > ^(3770) — > D + D~ , close to the threshold, and are thus almost at rest. This 
difference of production environments is important for observation of events from the decay 
D + — > K s tt + , which has a large rate and contributes to the same final state. These events 
are easily removed in the fixed target experiments by requiring all three charged pions to be 
consistent with a common vertex, and its residual contribution was estimated to be small. 
We are forced to take a different approach as the lower momentum K s does not produce 
clearly detached vertexes when K s — > tt + tt~ . Nevertheless we are able to clearly isolate the 
K s tc + channel, using the tt + tt~ invariant mass. 

Our analysis compares several different models for this decay, attempting to find the best 
description. One is an isobar model where we have included the best description of the a 
from Ref. jg] and the Flatte parameterization for the threshold effects on the /o(980) [8|. We 
use two other S wave models, both of which satisfy chiral constraints and respect unitarity. A 
model by Schechter and his collaborators (Schechter model) [sj is based on the linear sigma 
model of the chiral symmetric Lagrangian. It includes only the lowest lying tttt S wave 
resonances, the a and the /o(980). A model by Achasov and his collaborators (Achasov 
model) [13] is field-theory based and has been developed to describe scattering experiments. 
We compare the results of these three models of the resonance contributions to the Dalitz 
plot to see if one description is superior to the others and to understand differences among 
the models. 
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TABLE I: A comparison of the observed fit fractions in % from previous studies of D + — ► ir~Tr + ir + . 
The sum of all fit fractions is not necessarily equal to 100% due to the ignored interference terms. 
The "S wave 7r + " entry for E791 is the sum of the three entries above it. 



Mode 


E791 [3] 


FOCUS [4] 


0"7T + 


46.3 ± 9.2 




/o(980)tt+ 


6.2 ± 1.4 




/o(1370)vr+ 


2.3 ± 1.7 




S wave 7r + 


54.8 ±9.5 


56.0 ±3.9 


P°{770)tt+ 


33.6 ±3.9 


30.8 ±3.9 


/ 2 (1270)vr+ 


19.4 ±2.5 


11.7 ± 1.9 



In Section [TT] we briefly describe the CLEO-c experiment and the basic algorithms of event 
reconstruction. In Section II I II we describe the event selection for the Dalitz plot analysis. 
The formalism of fitting the observed Dalitz plot, and systematic cross-checks are given in 
Section [TV] Appendix IVHI describes in detail the two 7r + 7r~ S wave models that we use, 
some of which are extensions of published theoretical work. We summarize our results in 
Section [V] 



II. DETECTOR AND EXPERIMENTAL TECHNIQUE 

CLEO-c is a general purpose detector which includes a tracking system for measuring 
momenta and specific ionization of charged particles, a Ring Imaging Cherenkov detector to 
aid particle identification, and a Csl calorimeter for detection of electromagnetic showers. 
These components are immersed in a magnetic field of 1 T, provided by a superconducting 
solenoid, and surrounded by a muon detector. The CLEO-c detector is described in detail 
elsewhere (ll| . 

This analysis utilizes 281 pb _1 of data collected on the ?/>(3770) resonance at 
■t/s ~3773 MeV at the Cornell Electron Storage Ring, corresponding to production of about 
0.78 x 10 6 D + D~ pairs. We reconstruct the D + — > 7r~7r + 7r + decay using three tracks 
measured in the tracking system. Charged tracks satisfy standard goodness of fit quality 



requirements |12|]. Pion candidates are required to have specific ionization, dE/dx, in the 
main drift chamber within four standard deviations of the expected value for a pion at the 
measured momentum. Tracks coming from the origin must have an impact parameter with 
respect to the beam spot (in the plane transverse to the beam direction) of less than 5 mm. 
We do not reconstruct the Kg — > 7r + 7r~ vertex, but the requirement on pion track impact 
parameter removes ~60% of events with Kg — > 7r + 7r~ decays. The remaining events from 
D + — > K s tt + represent about one third of those selected for the Dalitz plot. 

III. EVENT SELECTION 

Selection of events from the D + — > it~-K + it + decay is done with two signal variables: 

AE = E D — -Ebeam, (1) 



m BC = V^beam ~ Pd, ( 2 ) 
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FIG. 1: The mec vs - AE distribution of 
events passing all selection requirements de- 
scribed in the text. The center box shows the 
signal region for the Dalitz plot analysis. The 
two hatched boxes show the sidebands. The 
vertical and horizontal lines restrict the re- 
gions of events plotted in Figs. [2] and EJ re- 
spectively. 



FIG. 2: The mec distribution of events from 
the \AE\ < 2a (AE) range. Dashed curve 
shows a contribution from the background, 
dotted curve is a Gaussian part of the Crystal 
Ball function for the signal shape, and solid 
curve is total, signal plus background. Events 
between arrows are selected for the Dalitz plot 
analysis. 



where -Ebeam is a beam energy, and E D and are the energy and momentum of the 
reconstructed D meson candidate, respectively. The beam crossing angle of ~4 mrad is 
used to calculate the D meson candidate energy and momentum in the ^(3770) center of 
mass system. We require |Ai?| < 2a (AE), \niBc — m o\ < 2cr(raBc)> where resolutions 
a(AE) = 5.5 ± 0.4 MeV and a(m BC ) = 1-38 ± 0.03 MeV/c 2 represent the widths of the 
signal peak in the 2D-distribution shown in Fig. [TJ and the projections, Fig. [2] and Fig. 
To determine the efficiency we use a GEANT-based Monte Carlo simulation where one of 
the charged D meson decays in a signal mode uniformly in phase space, while the other 
decays to all known modes with relevant branching fractions. Simulated events are required 
to pass the same selection requirements as data. The shape of the background contribution 
in the Dalitz analysis is estimated using events from the two hatched side-band boxes shown 
in Fig. [TJ The sideband boxes are shifted in AE to select the background events whose 
7r~7r + 7r + invariant mass range is consistent with the signal box. 

This selection gives 6991 events in the signal box. From a fit to the ui-qq distribution, 
shown in Fig. [2J we find 2159±18 of these to be background. The K s — > 7r + 7r~ contribution 
to the sample of events in the signal box is easily seen as a sharp peak in the invariant 7r + 7r~ 
mass spectrum shown in Fig. [H The Kg contribution is well described by a Gaussian shape 
with resolution a(m 7T + 7T ^) = 3.5 MeV/c 2 both in data and the simulation. The number of 
events in the K s peak is 2239±77 from a fit to a Gaussian signal plus linear background. 
Excluding K s tt + fraction and the background leaves ~2600 signal events of the D + — > 
7r~7r + 7r + decay. From these yields we calculate branching fractions, B(D + — > 7t~7t + tt + ) = 
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Data: D+ -> n + n + n~ 



Data: D + -> n + n + n ' 



2000 





A E (GeV) 



0.45 0.5 0.55 

m(7r + n -) Low (GeV/c2) 



FIG. 3: The AE distribution of events from 
the \m-QC ~ m D\ < 2<r(mBc) range. Events 
between the arrows are selected for the Dalitz 
plot. 



FIG. 4: The m(7r + 7r~)Low distribution of 
events pre-selected for the Dalitz plot. A 
clear signal for the Kg — > ir + it~ decay is ob- 
served. Events in the range between the ar- 
rows, 0.2 < m 2 (7r+7r-) Low < 0.3 (GeV/c 2 ) 2 , 
are discarded from the Dalitz plot analysis. 
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FIG. 5: The Dalitz plot for D+ 
candidates. 



7T 7T + 7T + 



FIG. 6: The adaptive binning scheme. 
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(0.33 ± 0.01)% and B(D + -> K s tt+) = (1.59 ± 0.06)% (statistical errors are shown only), 
which are consistent with recently published CLEO-c results B(D + — > 7r~7r + 7r + ) = (0.34 ± 
0.02)% and B(D + -> K° s n + ) = (1.55±0.05±0.06)% Q. This cross-check demonstrates 
the quality of our simulation and validity of assumptions about the background level. 

The presence of two 7r + mesons impose a Bose-symmetry of the tt~7t + tt + final state. The 
Bose-symmetry when interchanging the two same sign charged pions is explicitly accounted 
for in our amplitude parameterization. We analyze events on the Dalitz plot by choosing 
x = m 2 (7r + 7r~)Low and y = wi 2 (7r + 7r~)Hi g h as the independent (x,y) variables. The third 
variable z = m 2 (7r + 7r + ) is dependent on x and y through the energy- momentum balance 
equation. This choice folds all the data into the top half of the kinematically allowed region, 
as is shown in Fig. [51 The contribution from D + — > KgTT + is clearly seen as the narrow 
vertical band with m(7r + 7r~)Low — %«• In our Dalitz plot analysis we do not consider 
events in the band 0.2 < m 2 (7r + 7r~) Low < 0.3 (GeV/c 2 ) 2 , which is approximately ten times 
our Kg — > 7r + 7r~ mass resolution. This leaves 4086 (signal and background) events for our 
Dalitz plot analysis. 



IV. DALITZ PLOT ANALYSIS 
A. Formalism 



This Dalitz plot analysis exploits the techniques and formalism described in Ref. [14] that 
have been applied in many other CLEO analyses. We use an unbinned maximum likelihood 
fit that minimizes the sum over N events: 

N 

£ = -2^1ogP(x n ,y n ), (3) 

n=l 

where V(x, y) is the probability density function (p.d.f.), depends on the event sample to be 
fit, 

{e(x, y) for efficiency; 

B(x,y) for background; (4) 

f sig Af s \M{x,y)\ 2 e{x,y) + (1 - f sig )Af B B(x,y) for signal. 

The shapes for the efficiency, e(x,y), and background, B(x,y), are explicitly x — y sym- 
metric, third order polynomial functions. To account for efficiency loss in the corners of 
the Dalitz plot, due to low momentum tracks that are not reconstructed, we use three 
multiplicative threshold functions that drop the efficiency to zero when one of the Dalitz 
variables x, y, or z is at their maximum values. The background shape parameterization 
also includes the non-coherent addition of three resonances p(770), ^(1270), and Kg. The 
signal p.d.f. is proportional to the efficiency-corrected matrix element squared, \Ai(x,y)\ 2 , 
whose fraction is / S i g . We estimate / S j g = 0.548 ± 0.013 from the fit to the mec mass 
spectrum after removing events of the K s contribution. The background term has a rela- 
tive (1 — / S ig) fraction. The signal and the background fractions are normalized separately, 
1/A/s — / \M.{x, y)\ 2 £(x, y)dxdy, I/Mb = / B(x,y)dxdy, which provides the overall p.d.f. 
normalization, / V(x,y)dxdy = 1. The matrix element is a sum of partial amplitudes, 

M=J2cbA r Q r F r , (5) 

R 
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where Ar is a mass and spin-dependent function, Qr is an angular distribution 14], and Fr 



is the Blatt-Weisskopf angular momentum barrier-penetration factor [15|. In our standard 
fit the complex factor cr = aRe l ^ R is represented by two real numbers, an amplitude clr and 
a phase <$r. These are included in the list of fit parameters and can be left to float freely or 
fixed. 

For well established resonances, such as p(770), / 2 (1270), / (1370), / (1500), / (1710), 
etc., Ar is modeled with the Breit-Wigner function 

1 

A R (m) = — : , (6) 

m R — wr — im R i R\m) 

where m is the 7r + 7r~ invariant mass, rriR and r^(m) are the resonance mass and mass 



dependent width [14j, respectively The Ar parameterization of the /o(980), whose mass, 
mj , is close to the KK production threshold, uses the Flatte jl] formula 

-4/o(980)M = — 2 — 2 J7~2 1 J—\ i ~2 1 TZTVI' ( 7 ) 



m 



fo -m 2 - i[g 2 fo7T7T p n7T {m) + ff* fokF p kkM] 

where gf nn and 9^kk are ^ ne /o(980) coupling constants of the resonance to the ttti and KK 
final states, and p a b{m) = 2p a /m is a phase space factor, calculated for the decay products 
momentum, p a , in the resonance rest frame. 

We model a low mass tctt S wave, a or /o(600), in a number of ways. To compare our 
results with E791 we try a simple spin-0 Breit-Wigner. We also tested a complex pole 
amplitude proposed in Ref. [6j: 

Mm) = 2 1 2 ' ( 8 ) 

where m a = (0.47 — z0.22) GeV is a pole position in the complex s = m 2 (7r + 7r~) plane 
estimated from the results of several experiments. We also consider two comprehensive 
parameterizations of the low mass tttt S wave. One of them, suggested by J. Schechter, is 
discussed in Section IIV C\ and its formalism is presented in Appendix IVII A[ Another one, 
suggested by N.N. Achasov, is discussed in Section HVDt and its formalism is presented in 
Appendix EHEJ 



B. Fits with Isobar Model 

We begin our Dalitz plot analysis by attempting to reproduce the fit results E791 jl]. 
Our amplitude normalization and sign conventions are different from E791. We therefore 
compare the phases and fit fractions only. In Fit^l the contributions from p(770)7r + , 
/o(980)7r + , / 2 (1270)7r + , /o(1370)7r + , p(1450)7r + , and non-resonant intermediate states are 
included. Fit#l gives a probability of ~ 0. We checked that the inclusion of a air contri- 
bution, Fit #2, agrees better with the data giving a fit probability of ~ 20%. We obtain 
good agreement comparing our results with Fit#l and Fit#2 discussed in Ref. p. Then, 
we systematically study possible contributions from all known 7r + 7r~ resonances listed in 



Ref. [16|]: p(770), / 2 (1270), / (1370), p(1450), / (1500), /o(1710), and / (1790). We do not 
consider /2(1525) due to its negligible branching fraction to 7r + 7r~. We assume that high 
mass resonances ps(1690) and p(1700), having non-uniform angular distributions at the 
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edge of the kinematically allowed region, are well enough represented by /o(1710), which is 
a KK dominated resonance. The asymptotic "tails" of other known higher mass resonances, 
/2(1950), /4(2050), are effectively accounted for in our fits by the /o(1790) contribution. We 
also include a unitary amplitude parametrization of the n + ir + S-wave with isospin 1=2 from 
Ref. (l7| . For the /o(980) we use the Flatte formula, Eq. [7J with parameters taken from the 
recent BES II measurement [llj]. For the a we switch to a complex pole amplitude, Eq. [HI 
rather than the spin-0 Breit-Wigner used by E791. 

Starting from the contributions clearly seen in our fit, which is equivalent to Fit#2 
of E791 [3j, we add or remove additional resonances one by one in order to improve the 
consistency between the model and data. We use Pearson's \ 2 statistic criterion [l6[ for 
adaptive bins to calculate the probability of consistency between the p.d.f. and the data 
on the Dalitz plot. The bins are shown in Fig. [6j We also consider the variation of the 
log likelihood to judge improvement. We keep a contribution for the next iteration if its 
amplitude is significant at more than three standard deviations and the phase uncertainty is 
less than 30°. Table HT1 shows the list of surviving contributions with their fitted amplitudes 



TABLE II: Results of the isobar model analysis of the D + — > 7r~7r + 7r + Dalitz plot. For each 
contribution the relative amplitude, phase, and fit fraction is given. The errors are statistical and 
systematic, respectively. 



Mode 


Amplitude (a.u.) 


Phase (°) 


Fit fraction (%) 


p(770)vr+ 
/o(980)vr+ 
/ 2 (1270)vr+ 
/ (1370)vr+ 
/o(1500)^+ 
a pole 


1 (fixed) 
1.4±0.2±0.2 
2.1±0.2±0.1 
1.3±0.4±0.2 
1.1±0.3±0.2 
3.7±0.3±0.2 


O(fixed) 
12±10±5 
-123±6±3 
-21±15±14 
-44±13±16 
-3±4±2 


20.0±2.3±0.9 
4.1±0.9±0.3 
18.2±2.6±0.7 
2.6±1.8±0.6 
3.4±1.0±0.8 
41.8±1.4±2.5 



and phases, and calculated fit fractions. The sum of all fit fractions is 90.1%, and the fit 
probability is ~28% for 90 degrees of freedom. The best p.d.f. and the two projections of 
the Dalitz plot and selected fit components are shown in Figs. [7J [BJ and[9j For contributions 
that are not significant we set upper limits at the 95% confidence level, as shown in Table IrLTl 
The "N.R." represents a non-resonant contribution which is assumed to populate the Dalitz 
plot uniformly with a constant phase. 

TABLE III: Upper limit on the fit fraction, at the 95% confidence level, for contributions that we 
do not find significant in the D + — > 7r~7r + 7r + isobar model Dalitz plot analysis. 



Mode 


Upper limit on fit fraction (%) 


p(1450)vr+ 


<2.4 


N.R. 


<3.5 


1=2 tt + it + S wave 


<3.7 


/o(1710)tt+ 


<1.6 


/ (1790)7T+ 


<2 



The systematic uncertainties, shown in Table [TT1 are estimated from numerous fit vari- 
ations. We study the stability of the nominal fit results by adding or removing degrees of 
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FIG. 7: The signal p.d.f. for the isobar model 
fit described in the text. 




FIG. 8: Projection of the Dalitz plot onto 
the m 2 (ir + 7r~) axis (two combinations per D + 
candidate) for CLEO-c data (points) and iso- 
bar model fit (histograms) showing the vari- 
ous components. 



FIG. 9: Projection of the Dalitz plot onto the 
m 2 (7r + 7r + ) axis for CLEO-c data (points) and 
isobar model fit (histograms) showing the var- 
ious components. 
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freedom, varying the list of contributions to the Dalitz plot, changing the event selection, 
and varying the efficiency and background parameterizations. The systematic uncertainty 
of each fit parameter is estimated as the quadratic sum of the mean and root mean square 
values of the distribution of the changes in the parameter from its value in the nominal 
fit. For example for the poorly established resonances / (980), / (1370), and a pole, we 
allow their parameters to float and the variations of the other fit parameters contribute to 
the systematic errors. The nominal and fitted values of these parameters are presented in 
Table IIVI The fit results when the parameters are allowed to float do not vary from the 
nominal values by more than two standard deviations. 

TABLE IV: Parameters for the poorly established resonances used in the nominal isobar model fit 
and their fitted values when they are allowed float. 





Parameter 


Nominal Value 


Fitted Value 


Signal fraction 


/ sig from Eq. H 


0.548 


0.552±0.020 


/o(980) 


"Vo(980) (MeV/c 2 ) 


965 


953±20 




g fo7Tn (MeV/c 2 ) 


406 


329±96 






2 - fixed 


2 - fixed 


/o(1370) 


«Vo(i370) (MeV/c 2 ) 


1350 


1259±55 




r/ (i370) (MeV/c 2 ) 


265 


298±21 


a pole 


Re{m a ) (MeV/c 2 ) 


470 


466±18 




Im(m a ) (MeV/c 2 ) 


-220 


-223±28 



C. Schechter Model 



The isobar model drawbacks are most apparent in the S wave tt + tt~ sector where wide 
resonances overlap and unitarity is not fulfilled. The model of Joseph Schechter and co- 
workers in Refs is based on the meson part of the chiral invariant linear sigma 



model [20j Lagrangian. Poles are handled using K-matrix regularization which respects 
unitarity by definition. Details of the parameterization are discussed in Appendix IVII A[ 
and here we only summarize the meaning of the fit parameters. 

In our isobar model Dalitz plot fit the 7r + 7r~ S wave is represented by a complex pole 
for the a, the Flatte for the / (980) and two Breit-Wigner for the / (1370) and / (1500). 
Schechter's S wave amplitude, Eq. [221 (Appendix IVII Aj) . parameterizes simultaneously the 
a mixed with the /o(980) in strong and weak interactions. The Schechter model describes 
the mixed a and /o(980) contributions to the Dalitz plot with seven parameters: the bare 
masses m a and m/ ; the strong mixing angle ip between the a and /o(980); the total S wave 
amplitude asw an d phase 4>sw] an d the relative weak amplitude aj and phase 4>f of the 
/o(980) with respect to the a amplitude. A combination of these parameters in the model 
gives the total 7r + 7r~ scattering phase, 8(m), and an overall S wave amplitude, A sw , for the a 
and / (980) contributions. Operationally we replace the isobar a and / (980) contributions 
by the function of Eq. [22] times csw — a>sw e w • The Breit-Wigner's parameterization is 
still used for the / (1370) and / (1500). 

In an initial fit #S1, shown in Table [V] we fix all amplitudes and phases to their values 
from our isobar model fit, fix the S wave model parameters as in Eq. [191, float the S wave 
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TABLE V: S wave amplitude parameters in the fit of the Schechter model described in the text to 
the D + — > 7r~7r + 7r + Dalitz plot. 



Mode 


#S1 


#S2 


#S3 


m a (MeV/c 2 ) 


847 


758±36 


745±55 


m /o (MeV/c 2 ) 


1300 


1385±101 


1221±128 


^ (°) 


48.6 


45±5 


38±9 


asw 


4.1±0.2 


3.9±0.4 


4.5±0.6 


<Psw (°) 


54±3 


54±4 


55±6 


°/o 


3.8±0.2 


4.2±1.5 


2.1±1.5 


</>/o (°) 


23±3 


22±5 


21±5 


FF(S wave) 


45.9±1.9 


46.4±4.8 


43±12 


Ei^i (%) 


92.1 


90.6 


88.3 


Pearson/iVd.o./. 


116.3/96 


100.4/93 


99.6/87 


Probability (%) 


7.8 


28.2 


16.8 


-2£logL 


414 


398 


397.3 



amplitude asw an d phase <psw, an d noa t the relative / (980) amplitude a,f and phase 0/ o 
in Eq. [22] This fit gives a probability of 8% which indicates the Schechter model for the S 
wave is an acceptable description of the data. 

In a second fit, #S2 in Table \V\ we start from the parameters obtained in #S1 and allow 
the bare masses m a , rrif , and the strong mixing phase ip in Eq. [22] to float. This fit gives 
a probability of 28% and m a = (758 ± 36) MeV/c 2 , which is ~ 3 standard deviations lower 
than the values obtained in Ref. as also shown in our Eq. [T9l The mass mj and the 
phase ip are statistically consistent with the results in Ref. [l9j |. 

Fits #S1 and #S2 are used for an initial assessment of the Schechter S wave parameters 
relative to the isobar model fit. In a final fit, #S3 in Table \V\ we float the Schechter S wave 
model parameters and all the parameters of the other contributions. The results of fit #S3 
are shown in Figs. [TU1 and [TT1 in projections of the Dalitz plot. Figure [121 shows the isolated 
S wave contribution to the Dalitz plot, and Fig. [TBI shows the tot scattering phase, 5(m), 
defined in Eq. [2D] in Appendix IVH Al The total signal contribution is very similar to that 
shown in Fig. [71 Figure dH shows the complex amplitude Asw from Eq. [22] as the real and 
imaginary parts, the magnitude and complex phase. 

Employing the Schechter model changes the fit parameters for the non S wave contribution 
by less than the systematic uncertainties in the isobar model fit. We also note that the 
amplitude and fractions of /o(1370) and /o(1500) tend to be larger in the Schechter model 
fit. This model gives an acceptable fit probability ~17% when it is used to describe the a 
and /o(980) fractions in our data. The S wave fit fraction, (43±12)%, is consistent with a 
sum of fit fractions from a, (41.8±1.4±2.5)%, and / (980), (4.1±0.9±0.3)% in the isobar 
model. We find the Schechter S wave model parameters, listed in Table [Vj are consistent 



with the values in Ref. [19j. Our data are consistent with both the isobar and Schechter 
models. 
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X+Y projections 
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FIG. 10: Projection of the Dalitz plot onto 
the m 2 (7r + 7r~) axis (two combinations per 
D + candidate) for CLEO-c data (points) and 
Schechter model fit #S3 (histograms) showing 
the various components. 



FIG. 11: Projection of the Dalitz plot onto the 
m 2 (-7r + 7r + ) axis for CLEO-c data (points) and 
Schechter model fit #S3 (histograms) showing 
the various components. 
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FIG. 12: The isolated S wave contribution of 
Schechter model fit #S3 on the Dalitz plot. 



FIG. 13: The tt + tt~ scattering phase S(m), 
Eq. [20l calculated for parameters from 
Schechter model fit #S3 to the D + — > 
7r~7r + 7r + Dalitz plot. 
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FIG. 14: Complex S wave amplitude from Schechter model fit #S3 to the D + — > 7r~7r + 7r + Dalitz 
plot. The real and imaginary parts, the magnitude and phase are shown as a function of 7r + 7r~ 
mass. 



D. Achasov Model 



In Refs. |2l[-|25j and references therein, a irir S wave interaction is studied for tt-k — > tt-k, 
lm — > KK, (f) — > (/o — (7)7 — > 7T7T7, and 77 — > 7r7r processes in a manner motivated by field 
theory. The tttt S wave production and the final state interaction (FSI) mechanism in D 
meson three-body decays have not yet been considered in the framework of this model. In 
Ref. [l(| the tttt S wave amplitude in D + — > ti^ti + ti + decay is discussed. The developed 
formalism is described in Appendix IVIIBt and here we only summarize the meaning of the 
fit parameters. The Achasov model treats the tttt S wave contribution to D + — > 7t~7t + tt + 
via the sum of a number of amplitudes. There is a contribution from the non- resonant, 
point-like 7r~7r + 7r + production amplitude; direct resonance production via the D + — > a7i + , 
D + — > / (980)7T + ; and the rescattering terms from several intermediate states, 7r + 7r~, 7r 7r°, 
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and KK, to the final 7r + 7r~ state. Our parameterization has an amplitude, a and 
phase, 4>d+r,it+, for the direct resonance production term, accounting for the a and f com- 
ponents controlled by the coupling constants go+a-K+ an d fi , D+/ 07 r+- The contributions from 
rescattering have amplitudes and phases parametrized by a mo d e and mo dc plus a parameter 
from loop diagram contributions, <i mo de- We explicitly fit for the "mode" = 7r + 7r~, 7r°7r°, 
and KK rescattering contributions. The contribution from non-resonant 7r~7r + 7r + is also 
accoonted for the relevant point-like production amplitude parameter. 

We start with the parameters, shown in Table [III where the a pole and /o(980) are 
replaced by the S wave amplitude from Eq. [751 We fix all resonance parameters from our 
isobar model fit and float different sets of S wave parameters to assess their range. In four 
fits we float the amplitude, a mo d e , phase, mo de, and the offset parameter, d mo d e , (or coupling 
constants gD+aw+ an d gD+f w+ in case of direct cr or f meson production) for sub-modes 
7r + 7r~ — > 7r + 7r~, KK — > 7r + 7r~, 7r°7r° — > 7r + 7r~, or "DRtt" , respectively. For each of the single 
sub- modes we get a fit inconsistent with data. In five fits we float a mo d e , mo de, <^mode (or 
gi)+aTv+ an d 5 , d+/ 7t+) parameters for each combination of two sub-modes. All fits without 
the 7r°7r° — > tt + tt~ sub-mode show probability of consistency with the data ~10%, while 
models with the 7r°7r° — > 7r + 7r~ sub-mode are poorly consistent with the data. In three fits 
we include three or more sub-modes. These have a consistency with the data of ~10%, but 
give poor statistical significance for the amplitude parameters. Fit $A1 allows full freedom 
for all the S wave sub-modes and gives a probability of consistency with the data of ~19%, 
with 2-3 standard deviation significance for the amplitude parameters. Its results are shown 
in Table [VD 

We begin again with parameters of Fit #A1 and float or set to zero amplitude the 
parameters of the / 2 (1270), / (1370), and / (1500) contributions from our isobar fit. In 
Fit #A2 we float all the S wave parameters and all resonance parameters for the i /2(1270), 
/o(1370), and /o(1500) contributions. Variations of the nominal fit parameters, shown in 
Table [TTl are within the range of the isobar model uncertainties. Fit $A3 is like Fit #A2, 
but the contributions from / (1370) and /o(1500) scalar resonances are set to zero. The fit 
quality change from Fit $A2 to Fit $A3 is small. The S wave of the Achasov model has 
enough freedom to substitute for the contribution of the /o(1370) and /o(1500) resonances. 
The results of these two fits are shown in Table IVIl The results of Fit $A2 are shown 
in Figs. [ToTfTTI giving the Dalitz plot projections onto the m 2 (7r + 7r~) and m 2 (7r + 7r + ) axes, 
and the representation of the S wave complex amplitude. Our data are consistent with the 
isobar, Schechter, and Achasov models. 

E. Discussion of Models 

We have tested three models of the low mass 7r + 7r~ S wave in D + — > 7r~7r + 7r + , and we 
find little variation of the parameters describing non S wave contributions. The fit gives 
similar S wave contributions for all three models. We show this by plotting the relevant 
complex functions describing the S wave. Figure [18] shows the Flatte and the complex- 
pole parameterizations for / (980) and a, respectively, for our isobar model fit to the data. 
Figure [TH shows the results of the Schechter model fit, and Fig. [T7J shows the results of 
the Achasov model fit. In Figs. [191 EQ] we compare the ttti S wave amplitude and phase 
in the accessible mass region from threshold to 1.7 GeV/c 2 for these three models. The 
solid curve corresponds to the Schechter model fit to our Dalitz plot, the dashed curve is for 
Achasov model fit, and the ±lcx of the amplitude and phase parameters range of the S wave 
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TABLE VI: Fit results for the Achasov model as described in the text. 



Sub- amplitude, 


#A1 


#A2 


#A3 


parameters 








DRir 








a D+Rn+ 


1-fixed 


1-fixed 


1-fixed 


0D+Rn+ (°) 


-3±32 


-66±7 


-92±13 


9D+U-K + 


24±11 


39±8 


21±12 


9D+f w+ 


27±11 


267±24 


132±44 


TT^TT^ — > 7T + 7r = ' = 








@>7T7V 


0.25±0.08 


0.31±0.04 


0.25±0.07 


<f>™ (°) 


104±12 


70±9 


93±9 




1.5±0.3 


2.2±0.2 


2.9±0.3 












0.56±0.39 


1.35±0.15 


1.80±0.4Q 


(°) 


110±24 


107±7 


81±12 




0.02±0.21 


0.90±0.09 


0.37±0.10 


7r°7T — > 7T + 7r~ 








a 7r°7r° 


0.13±0.07 


0.11±0.03 


0.06±0.05 


07r%°(°) 


41±31 


149±23 


0±41 






^7T7r 




Fit fractions (%) 








EiFFi 


112.3 


140.4 


117.1 


7T + 7r~, 2X 


32.1±9.8 


37.5±3.6 


34.2±5.3 


7T + 7T + 


6.1±5.0 


16.6±3.2 


9.9±3.0 


Fit goodness 








Pearson/rid.oj. 


100.7/89 


96.9/83 


106.8/87 


Probability (%) 


18.7 


14.1 


7.3 


-2£logL 


398.6 


394.7 


405.1 



contribution in the isobar model is indicated by the two dotted curves. The S wave shapes 
are quite similar up to the interplay with other resonances, and with the data set we have 
in hand we are not sensitive to the details of the S wave parameterization. 

V. SUMMARY 

Using a sample of 0.78 million e + e~ — > ^(3770) — > D + D~ events collected in the CLEO-c 
experiment, we performed a Dalitz plot analysis of the D + — > 7r~7r + 7r + decay. Our nominal 
results, obtained within the framework of the isobar model and shown in Table [Til reinforce 
the previous conclusion 0] , [ij that a sizable an + component is required, in addition to other 
intermediate states p(770)tt+, / 2 (1270)tt + , / (1370)tt + , / (1500)7r+, and /o(980>+, in order 
to describe the D + — > 7r~7r + 7r + decay. The systematic uncertainties are estimated by varying 
the fit parameters from their nominal values. We also show in Table IIVI a set of optimal 
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FIG. 15: Projection of the Dalitz plot onto 
the m 2 (7r + 7r~) axis (two combinations per 
D + candidate) for CLEO-c data (points) and 
Achasov model fit #A2 (histograms) showing 
the various components. 



FIG. 16: Projection of the Dalitz plot onto the 
m 2 {'K + iT + ) axis for CLEO-c data (points) and 
Achasov model fit #A2 (histograms) showing 
the various components. 



parameters for the a, / (980), and / (1370) resonances based on our isobar model fit to the 
D + — > 7r~7r + 7r + Dalitz plot. Limits on contributions from p(1450)7r + , non-resonant, 1=2 
7T + 7r + S wave, /o(1710)7r + , and / (1790)7r + , shown in Table HTI| are set at 95% confidence 
level. 

We tested other models of the low mass tttt S wave contributions and in each case ob- 
tain optimal parameters. In Table V] we summarize results for the model suggested by 
J. Schechter and co-workers model show consistent values for 

the parameters. We also apply the S wave model suggested by N.N. Achasov et al. fioj . 
This model has more freedom in sub-modes than we are confidently able to define with our 
data. Possible solutions are presented in Table IVII Further progress with this model can be 
achieved if several D meson decay modes with higher statistics are analyzed simultaneously. 

For all 7T7T S wave models we find that their fit fraction exceeds 50%, and confirm results 
of previous experiments of a significant contribution from a low mass 7t + tt~ S wave in the 
D + — > 7r~7r + 7r + decay. Table IVLT1 compares the fit fractions from the fits to the three models 
described above. The S wave fit fraction in Achasov model is three standard deviation larger 
than in the Isobar and Schechter model. The sum of all fit fractions is also larger in Achasov 
model, that indicates on difference in interference terms. The fit fractions for sub-modes 
are consistent between these three models. Figures [TH] and [2D] compare the amplitude and 
phase, respectively, for the 7r + 7r~ S wave contribution we have found in the three considered 
models. With our given data sample all three S wave parameterizations adequately describe 
the D + — > 7r~7r + 7r + Dalitz plot. 
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FIG. 17: Complex S wave amplitude from Achasov model fit #A2 to the D + — ► 7r~7r + 7r + Dalitz 
plot. The real and imaginary parts, the magnitude and phase are shown as a function of 7r + 7r~ 
mass. 
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FIG. 18: Complex S wave amplitude (complex pole for a and Flatte for /o(980)) from isobar model 
fit to the D + —* 7r~7r + 7r + Dalitz plot. The real and imaginary parts, the magnitude and phase are 
shown as a function of 7r + 7r~ mass. 



VII. APPENDIX: ALTERNATIVE MODELS OF THE vrvr S WAVE 
A. Formalism of the tt-k S wave suggested by J. Schechter 

A tree level tttt — > tttt scattering amplitude for two resonances a and a strongly-mixed 
with phase if> is given in Eq. 3.2 of Ref. 



t; 





tree 



cos 2 tjj 



a(s) + 



13(8) 



rat — s. 



+ sin 2 if) 



a[s) + 



m% — s_ 



(9) 
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FIG. 19: The tt + tt S wave absolute ampli- 
tude for different models. 



3540307-019 



TD 2 
CO 

CD 
CO 
CO 
.C 



- 












/ 

/ 

/ 







y 








I< 

c 

fi 

. . , . 


sobar model range 

.chechter S wave 

ichasov S wave 
. , . . i . 



0.5 



1.5 



Mass (GeV/c 2 ) 



FIG. 20: The it + it S wave phase for different 
models. 



TABLE VII: A comparison of the observed fit fractions (FF) in % in the three models of D + — > 
7r~7r + 7r + . For the "Isobar" column, the "Low S wave 7r + " entry is the sum of the two entries above. 



Mode 


Isobar 


Schechter #S3 


Achasov #A2 


<T7T + 


41.8 ±2.9 






/ (980)vr+ 


4.1 ±0.9 






Low S wave ir + 


45.9 ±3.0 


43.4 ± 11.8 


75 ±7 


/o(1370)tt+ 


2.6 ± 1.9 


2.6 ± 1.7 


3.2 ±0.7 


/o(1500)tt+ 


3.4 ± 1.3 


4.3 ±2.4 


4.0 ±0.8 


/ 9°(770)vr + 


20.0 ±2.5 


19.6 ± 7.4 


18.4 ±4.0 


/ 2 (1270)^+ 


18.2 ±2.7 


18.4 ± 7.4 


23.2 ± 5.0 


1=2 7r + 7r + S wave 






16.6 ±3.2 


Ei FF i 


90.1 


88.3 


140.4 



where 



a(s) 



1 - 



4mJ 



I 1 2 s 

\m„ — m„ 



5±2 



s — Ami 



m„ + s — Ami 



mi 



(10) 



3\/l- 



16ttF2 



(11) 



s is the 7r + 7r~ invariant mass squared, and ^=0.131 GeV are the pion mass and the 
decay constant, m a and m„ are the bare masses of two scalar resonances, and ip is a strong 
mixing angle. We use the original notation of Ref. 19J, the tilde is used for all parameters 
relating to the second scalar resonance, a, which in our case is associated with / (980). 
Equation [9] can be re- written as 

A + (12) 



rpO 

1 tree 
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where 

A = a cos 2 ip + a sin 2 ip, (13) 

B = (3- {m\ - s) cos 2 %l) + (3 ■ (m 2 - s) sin 2 ijj, (14) 

P=(ml-s){ml-s). (15) 

According to the Dyson equation for the nir scattering, Eq. 3.3 from Ref. 19J gives an 
expression for a total scattering amplitude through the tree amplitude: 

T°( s ) = T °° free( ^ . (16) 

1 — «Tq tree( s ) 

The scattering amplitude is a complex number, Tq (s) = |Tq (s)|e i<5 ^- ) , then the tree amplitude 
can be associated with the tangent of the scattering phase, 

T V e ( S )=tan5( S ), (17) 

and we get an expression for cos 5: 

1 P 
cos S = ^= — = -= 18 

VI + tan 2 5 ^jpi + (p • A + Bf 

Expression 5(s) = arctan (T® tree (sfj defines a scattering phase in the range [— f,f]- This 
phase S(s) has two discontinuities at s — m 2 and s = m 2 ^ for parameters taken from Ref. 19 



m a = 0.847 GeV/c 2 , m fo = 1.30 GeV/c 2 , ^ = 48.6°. (19) 
In order to remove discontinuities we add a phase-shift +ir above each bare mass: 

S(s) = arctan (T ° iree ( S )) + tt(#( S - m 2 ) + 9(s - m%)), (20) 



where 9(x) is a step function, that makes the phase smooth, as shown in Fig. 

In this model the production amplitude is obtained from the total scattering amplitude, 
Eq. HSl by replacing the first tree level ttty — > nn scattering diagram amplitude, T° tree , by 
the resonance propagator (m 2 — s) _1 with the coupling constant g a7T7T and keeping the proper 
re-scattering amplitude, represented by the "bubble sum" factor (1 — iT$ tree )~ l : 

^ = -^• 1 -to COB 6. (21) 

m 2 a - s 1 - tree m%-8 

Extending Eq. [21] (Eq. 15 from Ref. |9|]) for the case of two resonances a and /o(980) we get 
the total production amplitude with relative weak interaction mixing factor a/ e l ^o 



A-sw — A a + Af = cos 5 ■ e 



m^ — s m 2 j Q — s 



(22) 



Note, that Eq. [22] does not contain singular terms because both poles are contracted into 
the P factor from cos S, Eq. [15] and [18] For the first iteration we set 

o A = l, 0/ o = O o . (23) 
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It should be noted, that in the frame of this model, a is a scalar tttt resonance which has 
a bare mass m a as a parameter. The bare mass does not coincide with a peak position as 
in case of Breit-Wigner, that is clearly seen in Eq. [19] for the mass of /o(980). This simple 
model does not take in to account that the scalar resonances may have other decay modes, 
coupled channels. For example, it is well known that /o(980) has a KK decay mode with 
a mass dependent rate as large as ~20%. Presumably, this amplitude, obtained from the 
chiral Lagrangian, works well in the region close to the production threshold. In the case of 
SU(3) symmetry it accounts for the two low mass resonances a and /o(980). Other higher 
mass resonances such as the /o(1370) and /o(1500) are not taken into account. These issues 
restrict the precision and limit the application of this model. 



B. Formalism of the tttt S wave suggested by N.N. Achasov 

1. D + — > tt~tt + tt + total amplitude 



In this section we summarize a suggested formalism [lOj for a parameterization of the tttt 
scalar amplitude in the D + — > tt~tt + tt + decay, and present the details of our implementation 
in the Dalitz plot fitter with some relevant cross-checks. For the D + — > tt~tt + tt + decay 
Ref. [13] suggests the use of a tttt S wave amplitude that is a superposition 



A{D+ - 


7T+7T+ TT ) = 


AP\D+ -> 






+B[D+ 


-> TT+(tT+ TT' 


-> TT+ TP") - 


» TT+TT+ 7T _ ] 


+ <-> 2] 


+E[D+ 


- (VT+7T+ - 




» 7T+7T+ 7T _ ] 




+F[D+ 




) ^TT+TT 


"] «-» 


2] 


+B[D+ 


-> tt+(tt tt - 


-> 


tt+tt+tt~] - 


<-> 2] 


+C[D+ 


-> tt+(K+K- 


- -> TT+ 7T-) 




] + C[l «-> 2] 


+C[D+ 


-^7T+(K a K° 


-> TT+TT-) - 


-> TT+TT+ 7T"] 


+ 2]. 



(24) 



of a point-like, A p ', direct resonance, F, and non-resonant production terms, B, C, E, 
followed by the re-scattering in to the tttt final state. Here we list the definitions of all the 
sub-amplitudes in Eq. |2"%1 

The point-like D + — > -KiTi^^ amplitude is associated with a constant a: 

A pl (D+ -> TTi+Tr+Tr-) = 167ra. (25) 

After the point-like production one would expect Tr^Tr - — > Tr^Tr - and ?T^7r _ — > tt^tt - scat- 
tering, which we parametrize as a mass dependent amplitude 

B[D + -> TT^ (tT^tT -> TT^TT") -> TT^ 7T^ 7T~, 771 = 77^-] 

'2„„, , 1. 



L w+w -(m|a,p) • (-T °(m) + -T 2 (m))- (26) 



Functions L^+ T -(m|a,p), Tq(t7i), and T 2 (m) are described below. An exotic 1=2 S wave 
7r+7r^ — > 7T+7r^ scattering is discussed in Ref. [lj 



.E[.D + — > (n^n^ — > TT^TrJ) — > Tr^Tr^Tr ,m = 771^+^+] = L 7r + 7r +(m|a, r) • TJ^m). (27) 
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It is assumed that the a and fo mesons can be produced directly in the D + — > tt + ct and 
D + — > 7r + /o decays (we use the u DRir" notation), with an amplitude of 

F[jD + ^ + /o) _> 7r + 7r + 7r - )m = m ^_] = T % e + V + (m). (28) 

The point-like D + — > 7r + 7r°7r° amplitude is associated with another constant a 

A pl (D + -> 7T + 7r°7r ) = 16™. (29) 

Subsequent 7r°7r° — > 7r + 7r~ rescattering may also contribute to the final state via the ampli- 
tude 

B[D + — > n^(n n — > 7r^7r~) — > Tr^TT^ - ,Tn = m^o^o] 

= L^o(m\a,q) ■ (Jr o °(m) - ~T 2 (m)Y (30) 

In the above equations we assume that q = r = p. 

The point-like production amplitudes for D + — > n + K + K~ and D + — > tt + K°K are 
represented by the constants c and c, 

A p '(£>+ -»■ n + K + K-) = IQttc. (31) 

A p '( J D + -> 7i + K°K°) = 16nc. (32) 
Then, two terms account for the relevant rescattering amplitudes K + K~ — > 7r + 7r _ and 

C[L> + -> u+{K + K- -> 7T+7T-) -> 7r5 l "7rj7r~,m = m^ +K -] 

= L^ +K -(m|c,s) -T°(K + K- -> 7r + 7r-,m), (33) 

+ - 7T+(K Z - 7T+7T-) - TT+TT+TT- m = m rf ] 

= L^ -o(m|c,t) • T °(K Z° - 7r + 7r-,m), (34) 

where we assume that offset parameters are equal, £ = s. 

In above equations we use the function L agi (m\c, d), which represents a contribution from 
the loop diagram 

L a a(m\c, d) = 16ttc • I W^™) + P<*( m H ln i+p°°M + d > m > 2m «» (35) 

I -|Poa(m)| + |poa("i)|f arctan |p oa (m)| + d, m < 2m a , 

where 

I Poo (HI = \l^m 2 a /m 2 - 1, (36) 

paa(H = ^/l - 4m 2 ./m 2 . (37) 

Below all definitions, required for parametrization of the amplitude in our case, are re- 
written from the recent Ref. \2l 
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2. T ° = r (7T7r -► 7T7T,m) 

Equation 23 from Ref. [25] gives the S wave amplitude Tq of 7r7r — > tttt scattering with 



o 

1=0 is 

T° , T °<™ ->7nr,ro) = = i-^ + . T » ~ (m) . (38 ) 



Equation 24 from Ref. [25j gives the total phase 

5° = <5 °(m) = S^(m) + 5 res (m). (39) 
Equation 25 from Ref. [25[ defines the resonant part of the S matrix 

S ° res (m) = r]° {m)e 2i5r ^ m) = 1 + 2ip mr (m) ■ T ° res {m), (40) 
which can be described by the inelasticity 

V ° (m) = \S° res (m)\ (41) 



and resonant phase 



i /%sr 



firesim) = - • arctan —) " ). (42) 

The chiral background shielding phase S]^(m), motivated by the a model, is taken as 
Eq. 26 from Ref. 

tan^ = -— (bo-h^ + b 2 ^) x - 1 N2/a2 , (43) 



V m*/ 1 + (2p 7T ) 2 /A 

where 2p n = ^m 2 — 4m^, and (1 + (2p n ) 2 / A 2 )^ 1 is a cutoff factor. The value of parameters 
bo, bi, 62, and A used in our fits are listed in Table IVnTl The background phase is derived 
from Eq. 1431 

8% = arctan [tan 8™]. (44) 



3. Resonance amplitudes Tq ^+ Rn + (ni) and Tq res (m) 

In Eq. [28] [38J, and HO] we use a brief notation for the production and scattering resonance 
amplitudes expressed through the mixing matrix operator G RR , 

rpO res _ JSV'(m) 9D+RTT+GRR'9R'ir+ir- , A r\ 

RR> 167F 

rpO res _ 9RmrG RR ,gR> nn 
RR> 

Note the difference between specific coupling constants and the exponential factor in Eqs. 05] 
and US 
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I T$(KK TTTT,m) 



The S wave amplitude of KK — ^ tttt scattering, taking in to account mixing through RR' 
resonances (i.e. a and /o(980) mesons) is given by Eq. 3 from Ref. [251 ]: 

9RK+K- G RR I gR'TV + TT- 



T°(K + K- 



7T 7T , m 



E 

RR' 

E 

RR' 



16tt 



9 rk qk (, Grr'9r'it+ 
16tt 



where Eq. 4 from Ref. [25J defines 



Equation 28 from Ref. [25 



is 

KK 



tan<5 B = f K (m) -2p K = f K (m) ■ ^Jm 2 - 4m 2 K , 
where Eq. 36 from Ref. [25j gives 



fici™) — ~ arctan 



m — m 



m\ 



A 



K, 



and we find the phase as 

5%x = arctan(tan5f^). 
The value of parameters mi, mi, and Ak used in our fits are listed in Table rVTTTl 



(47) 
(48) 

(49) 
(50) 

(51) 
(52) 



5. An exotic 1=2 amplitude Tq (m) = Tq (tt + tt + — > 7r + 7r H 



,m) 



According to Ref. [17| the 1=2 7r + 7r + — » 7r + 7r + the rescattering amplitude is given in a 
unitarian form 

mil \ — mil + + + + \ 'ov" e 1 /ro\ 
T (m) = T (7T + 7T + -> 7T + 7T + ,m) = . (53) 

The phase shift 5o( m ) * s parameterized by 



—a\/m 2 /A — m 2 
5 2 ( m ) = - V 7 d * . (54) 
ov 7 1 + bm 2 + cm 4 + dm 6 v ; 



From fit in Ref. [l7j to data for the 7r p — ^ n ir n process in Refs [27J and 28[], the parameters 
of Eq. [51 are a = (55.21 ± 3.18) deg/GeV, b = (0.853 ± 0.254) GeV" 2 , c = (-0.959 ± 
0.247) GeV" 4 , and d = (0.314 ± 0.070) GeV" 6 . 

The ?7o( m ) * s an inelasticity for the wave with total spin and isospin 2. In the mass 
range of m < m(pp) (~1.54 GeV) the inelasticity parameter rj^im) should be represente d by 
the smooth real function of m. An appropriate fit to data has been considered in Ref. [29], 
see their Fig. 2, and we use the approximation 

f 1, m < 1 GeV/c 2 

?7o( m ) = S oc cos— like smooth transition, 1 < m < 1.7 GeV/c 2 (55) 
[ 0.4, m > 1.7 GeV/c 2 . 

In our case we neglect the small D wave scattering amplitude Tf (7r + 7r + — > 
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6. Mixing matrix Grri (m) 



The mixing operator Grr> (m) is a matrix of inverse propagators, with rank equal to the 
number of mixed resonances. In case of mixing of two resonances R and R' this matrix has 
the form, following Eq. 5 of Ref. [25| . 



Grr> (m) 



D R (m) -Urr, (m) 
-U R/R (m) D R f(m) 



In general, the diagonal elements of this matrix are the inverse propagators 

Dr{tti) = m 2 R — m 2 — imT R (m), 



(56) 



(57) 



while the non-diagonal elements are polarization operators describing mixing. An expression 
for the inverse propagator of the scalar resonance is given in Eq. 6 from Ref. 25j , 



D R (m) 



m R — m 



Y,9Ra b [ReP a R\mR)-P a R \m)l 



(58) 



ab 



where J2 a b 9Rab[RePR{rn R ) — P R (m)] = Re[H R (m R )] — TlR(m) takes in to account the finite 
width correction. After Eq. 5 in Ref. [25] the non-diagonal terms of the polarization operator 
are given by equation 

iWm) = Y,9R>abP a R b {m) + CW, (59) 

ab 

where the constants C RR i take into account effectively the contribution of VV, 4P and other 
intermediate states and incorporates the subtraction constants for the R — > (PP) — > R' 



transitions. Here we use the notation from different publications, [221] — [25 



Pr(™) 



and 



16tt 2 K h RK ' 



UR{m)=^ a R b (rn). 

ab 



9R*tpab,, 

16tt 2 V 7 ' 



(60) 



(61) 



Eqs. 7-9 from Ref. [25[| (also Ref. [231, Eq. 30 and Ref. [24J, Eqs. 16,19,22) for m a > m b , 
m + = m a + m b , and m_ = m a — m b give 



P a \m) 



n In— +< 



Pab{m) 



m + In 



-7r\p ab (m) \ + 2\p ab (m) \ arctan 



-Pab{m) ■ In 



Pab(m) 



1 - 



1 - 



m: 



77T 



The constants g# a ft are related to the width, Eq. 11 from Ref. [25 

T(R^ab,m) = -^p ab (m). 



m > m + 

m_ < m < m + (62) 
m < m_ 

(63) 
(64) 
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7. Model parameters 



In the mixing operator Eq. [56] we account for seven intermediate states: tt + tt , 7r 7r°, 
K + K~ , K°~K°, rjr], rj'rj, and rj'r]'. We follow the conventions of Ref. [25} for coupling con- 
stants, motivated by the four-quark model. For the /o(980) and similarly for the o we 
use 

9f K°K° = 9foK+K-, gRn°n = 9Rtt+7t- /V2, gR ww = \J 3 /2g RlT + 7T - . (65) 

For the /o(980) coupling constants to rf^ff^ we use 



9fovv ~ 9fon'v' 



2^2 



-9f K+K~ , 9fov'v 



V2 



-9f K+K~ 



(66) 



For the a coupling constants to rf'^rft' we use 



9 



urfq 9<ril'r) 



-9, 



(TTT IT •> 



9 



3V2 



9c 



(67) 



Further we use the values of the parameters shown in Table IVHH which are taken from 
Fit 1 of Ref. [H. 



TABLE VIII: Achasov model parameters from Fit 1 of Ref. 25f] used in our calculations 



Parameter 


Value in Fit 1 [25] 


Parameter 


Value in Fit 1 [25] 


m fo , MeV 


984.1 




4.9 


m a , MeV 


461.9 


h 


1.1 


QfoK+K-, GeV 


4.3 


b 2 


1.36 


S/ott+tt-, GeV 


-1.8 


A, MeV 


172.2 


9aK+K~ i GeV 


0.55 


mi, MeV 


765.4 


SW+7r-> GeV 


2.4 


m 2 , MeV 


368.9 




-0.047 [26] 


A K , GeV 


1.24 



Check for 5 res (m) , S^(m), 5o(m), ^(m), etc. 



In order to check that the code for this parameterization works properly we reproduce 



plots from Ref. [25 



<5 res (m): We define the 5 res (m) as the phase of the complex function res (m) in Eq. H2J How- 
ever, this phase has discontinuities in the vicinity of each resonance mass, but not 
exactly at the resonance mass value. In further calculations we require that the phase 
is continuous, as shown in Fig. [22], by adding a phase shift of ti above each discontinuity 
point. This plot is consistent with Fig. 3 in Ref. [25| . 



^""(m): The background phase S^(m) is derived from Eq. HH as shown in Fig. [5TJ This plot 
is consistent with Fig. 2 in Ref. [25| . 
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FIG. 21: The background phase in tttt scat- 
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FIG. 23: The total -kit scattering phase, 
<5§(m), from Eq. ES 



FIG. 22: The phase of the resonance irir scat- 
tering, <5 res (m), from Eq. | 
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FIG. 24: The inelasticity, t]q, from Eq. I41I for 
277-^ < m < 2 GeV/c 2 . 



c5°(m): The total phase fflm) represented by Eq. EHis shown in Fig. [23 This plot is consistent 
with Fig. 4 in Ref . Q . 

T]q(jxi): The ^(m) derived from Eq. is displayed in Fig. I2~4l which shows that t]^(m) = 1 
at m < m K K confirming unitarity in tttt — > tttt scattering, consistent with Fig. 6 from 
Ref. @. 

We also tested all complex functions and their components from Eq. [2H In particular, 
Fig. [25] shows 5f R from Eq. [52j Fig. M shows 5 B from Eq. \M Figs. [271 ES show the loop 
integrals L K + K -(m\l, 0)/167r and L 7r + n -(m\l,0)/16Tr, respectively, from Eq. 
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FIG. 25: The background phase in KK scat- 
tering, 6g K , from Eq. [52] (solid curve), and 
its approximation by the phase space factor 
(dashed curve). 
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FIG. 26: The total background phase in 
KK — > -kit scattering, 5b = 5 B n + 5^ , from 
Eq.g9l 









3540307-027 














/ 
1 
1 

1 / 
\ f 


Im 








Re 






iii 




ii 


II 



0.5 1 1.5 2 



Mass (GeV/c 2 ) 



FIG. 27: The loop integral, 

L K + K -(m, 1, 0)/167T, from Eq. [35j The 
real (solid curve) and imaginary (dashed 
curve) parts of the complex function are 
shown. 



FIG. 28: The loop integral, 

L^+ w - (m, 1, 0)/167r, from Eq. [35j The 
real (solid curve) and imaginary (dashed 
curve) parts of the complex function are 
shown. 
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9. S wave implementation in the code of the Dalitz plot fitter 



As usually in a Dalitz plot analysis each amplitude fraction is taken with its own complex 
coefficient c mo d e = a mo d e e^ modo represented by two real numbers, an amplitude a mo d e and 
phase mo de- The loop integral in Eq. [33 has an additional offset constant <i mo de- Unitarity 
requires that (i mo de is real. All these constants, as well as unknown coupling constants 
g D+a7T+ and (?d+/ 7t+ from Eq. S3, are the fit parameters, which can be free to float or 
fixed. The actual parameterization for m = m x = min[m(7ri~7r~), m(7T2~7r~)] or m = m y = 
max[m(7r 1 ~7r _ ), m(7r^7r - )] is given by the amplitude 

A K +^-(m) = Wnc^ (68) 
+ L^ +7r -(m\c^,d^) ■ (-T °(m) + ^T 2 (m)^) (69) 

+ L n o n o{m\c n o n o,d T o T o) ■ f-T °(m) - ^T 2 (m)^) (70) 

+ L K+K -(m\c K +K-,d K + K -) ■ T%(K + K~ -> n + n~,m) (71) 

+ L K0 ^{m\c KOl?0) d KOl?0 ) ■ T°(K°K° -> tt+tt", m) (72) 

+ c D+R7r+ -T° r D ^(m). (73) 

The 1=2 7r + 7r + — > 7r + 7r + scattering amplitude for m = m z = m^fff^) is given by 

A 7r +^+(m) = L^+(m\c n7n d^) ■ T 2 (m). (74) 

It is worth noting that three terms in Eqs. [68], [69] and [71] have a common complex coefficient 
c n7T , appearing from the point-like term, and two of them have a common offset parameter 
d nn from the loop integral. The total contribution of Achasov's S wave in the Dalitz plot 
amplitude is 

Asw{m x ,m y ,m z ) = A w+7r -(m x ) +A 7r+n -(m y ) +A n + ir +(m z ). (75) 

The u DRn" sub-mode in Eq. [73] has a redundant freedom for amplitude factors due to the 
products CLD+R-K+ -gD+(j-K+ and a D +R 7T + •fi , D+/ 7r+- In our fits we fix a D + R7T + = 1, or a D + RlT + = 
to turn it off, and use coupling constants gD+a-w+ an d 5 , d+/ott+- 

For a first approximation we try to eliminate the number of free parameters in the func- 
tion. We assume d^o^o = d n7T and d R0 -^o = d^+x- from isospin symmetry. We note that the 

parameterization for K°K° — > tt + tt~ in Eq. [72]is nearly the same as that for K + K~ — > ix + ix~ 
in Eq. [7TJ The small difference appears due to the different masses of the K + and K° mesons. 
Keeping in mind this small difference between amplitudes we do not consider separate con- 
tributions from K°K — > 7r + 7r~ in this analysis. This means that the amplitude factor a K j^ 
includes both contributions from K + K~ — > tt + tt~ and K°K — > tt + tt~ . 

The amplitude for 7r°7r° — > 7r + 7r~ in Eq. [70] has a different isospin factor at TJ 2 compared 
to the amplitude for 7r + 7r~ — > 7r + 7r _ in Eq. [63 and different masses for tt° and tt + . In our fits 
we assume the equity d n o n o = d nn . The constant also accounts for the point-like term in 
Eq. [68] and is involved in 1=2 term, Eq. [74] that makes it different from the c^o^o. For this 
reason we consider the 7r°7r° — > 7r + 7r~ sub- mode separately from 7r + 7r~ — > 7r + 7r~. 
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